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The high-symmetry pseudocubic init cell is often used for modelling the electronic structure of halide per¬ 
ovskites. However, direct comparison of the band structure with more realistic low-symmetry phases is impeded 
by the zone folding. We utilize a spectral density technique in order to reduce the supercell band structure to a 
common Bloch basis. This allows us to compare the electronic structure of high- and low-symmetry phases as 
well as investigate effects of structural and compositional disorder on states near to the band edges that govern 
transport, dissociation and recombination of optical excitations. 
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I. INTRODUCTION 

A new class of organic-inorganic hybrid perovskite materi¬ 
als, such as (CH 3 NH 3 )Pbl 3 , has recently emerged as an alter¬ 
native choice for photovoltaic devices and immediately drew 
enormous attention M. (CH 3 NH 3 )Pbl 3 possesses a com¬ 
bination of characteristics that make it attractive for solar cell 
applications including a favourable band gap of 1.6 eV, ef¬ 
ficient optical adsorption, high level of mobility for charge 
carries of both polarities and exceptionally long carrier life¬ 
time i]. Device characteristics, such as the power conver¬ 
sion efficiency, can be further enhanced by using mixed halide 
perovskites, e.g., (CH 3 NH 3 )Pb(Cla^Ii_a ^)3 The demon¬ 
strated energy conversion efficiency of nearly 20 % com¬ 

bined with low-temperature solution processing technology 
hold a promise to lead the way towards low-cost solar energy. 

The crystal structure of (CH 3 NH 3 )Pbl 3 exhibits several 
polymorphs. At low temperatures, an orthorhombic 7 -phase 
prevails up to 161 K IH]. A tetragonal /7 structure is sta¬ 
ble in the intermediate temperature range of 161-330 K. The 
high-temperature cubic a-phase exists above 330 K |@]. Elec¬ 
tronic structure calculations performed for the cubic phase 
reveal a direct band gap of approximately 1.3 eV at the 
R(l/ 2 ,l/ 2 ,l/ 2 ) point of the Brillouin zone 1I9I-ITTI]. The di¬ 
rect gap moves to r(0, 0, O^s the cubic structure transforms 
into the tetragonal phase Cl [a, which has a more com¬ 
plex unit cell with larger volume. The transformation is ac¬ 
companied by widening of the band gap, which reaches the 
value of 1.6 eV as determined by photoluminescence mea¬ 
surements Cl and confirmed by ab initio calculations using 
hybrid functionals and quasiparticle GW approaches 

Further analysis of the band structure evolution in the 
course of the transition from high- to low-symmetry struc¬ 
ture is hindered by the Brillouin zone folding resulted from 
the increased size of the unit cell for y^-phase. Here we 
show that a direct comparison of the electronic structure be¬ 
tween a and f3 phases of (CH 3 NH 3 )Pbl 3 can be achieved 
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by projecting the band structure of /3-phase onto the Bril¬ 
louin zone of the cubic structure, which represents an ideal¬ 
ized version of low-symmetry phases. This analysis can fa¬ 
cilitate answering the following questions: (i) to which ex¬ 
tent does the Bloch character of band edges in the tetragonal 
phase inherit the R-character of the cubic phase, and (ii) what 
is the effect of dynamic disorder in the orientation of organic 
groups or substitutional alloying, e.g., in the mixed halide 
(CH 3 NH 3 )Pb(Cla;Ii_a ^)3 perovskite, on the Bloch states at the 
band edges? An uncertainty in the Bloch character of states 
near to the band edges in semiconductors can be used to reveal 
deterioration of charge transport characteristics due to carrier 
localization liMl, which can adversely affects the perfor¬ 
mance of solar cells. 


II. COMPUTATIONAL DETAILS 

The first-principles calculations were carried out using den- 
sity functional theory (DFT) in. ABINIT packageJmH^] 
and projector augmented-wave method were used 11 ^ with 
pseudopotentials adopted from a JTH library (version 0.2) 
. Perdew et al. construction of the generalized gra¬ 
dient approximation (GGA-PBE) was employed for the ex¬ 
change correlation functional. Completeness of the plane 
wave basis set was controlled though the cut-off energy, which 
was set at 15 Ha (408 eV) as determined from convergence 
studies. 

The lattice constants of (CH 3 NH 3 )Pbl 3 in pseudocubic 
phase (Fig. (Ha)) was optimized self-consistently by simulta¬ 
neous minimization of stresses and Hellmann-Feynman forces 
acting on atoms below 1 mHa/Bohr. The Brillouin zone was 
sampled using 8 x 8 x 8 Monkhorst and Pack ll24ll k-mesh. The 
calculated lattice constant of 6.47 A is slightly overestimated 
as compared to the experimental value of 6.28 A |@]. This 
result is consistent with typical trends for modelling of solids 
with GGA-PBE exchange-correlation energy functional iH. 

The structure of tetragonal phase was constructed by re¬ 
peating the cubic cell with the multiplicity of 2 x 2 x 2 
(Fig. [2a)). The octahedral tilting and orientation of methy- 
lammonium molecules were initially set based on previously 
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FIG. 1. Pseudocubic structure of (CH 3 NH 3 )Pbl 3 (a) and the corre¬ 
sponding relativistic band structure (b) along high-symmetry path in 
the Brillouin zone shown in the inset. 


reported structures of ^-(CH.sNH.s)Pbl8 resolved experimen¬ 
tally and theoretically 12^ [^. Then the supercell was fully 
relaxed. Calculated lattice parameters for tetragonal phase 
ao = 9.02 A and cq = 13.01 A are within 4% deviation 
from the corresponding experimental values at room temper¬ 
ature ao = 8.85 A and cq = 12.52 A 1^. It should be men¬ 
tioned that our tetragonal supercell differs from a conventional 
tetragonal unit cell by 45° rotation. The construction of the 
supercell by translation of the parent (cubic) structure along 
its lattice vectors allows us later to unfold the band structure 
and present the electron wave vectors in terms of the cubic 
reciprocal lattice. 

Unfolding of the band structure of supercells was per¬ 
formed by calculating the Bloch spectral weight using 
f old2Bloch packa ge [fo implemented in ABINIT. 

VE S T A 3 package was used for atomic structure visu¬ 

alization. 


III. RESULTS AND DISCUSSION 

We begin the analysis of the electronic structure of 
(CH 3 NH 3 )Pbl 3 with the high-symmetry pseudocubic struc¬ 
ture. The calculated relativistic band structure is shown in 
Fig-fflb). The band gap is direct near R(l/2,1/2,1/2) point 
of the Brillouin zone. The magnitude of the band gap is un¬ 
derestimated by approximately 1 eV in comparison to the ex¬ 


perimental value of 1.6 eV |[T3] . which is a typical result for 
DFT calculations of (CH 3 NH 3 )Pbl 3 when the bear GGA-PBE 
exchange-correlation functional is combined with spin-orbit 
coupling. The agreement can be further improved by using a 
modified Becke-Johnson potential 12^ . however it is not re¬ 
quired here due to the comparative nature of our study. A 
Rashba-like breaking of the spin degeneracy El is observed 
at the band ec^s (particularly strong in the conduction band). 
Zheng et al. 11^ reported a different chirality of the spin tex¬ 
ture at the top valence band and the bottom conduction bands 
that results in the long carrier lifetime due to constrains im¬ 
posed by the spin section rules. 

In order to elucidate origin of the Rashba splitting, the 
band structure calculation was performed for CsPbG in the 
pseudocubic structure (not shown here). Atomic position 
and lattice perimeters were kept at the equilibrium values of 
(CH 3 NH 3 )Pbl 3 , except for the methylammonium that was 
substituted by Cs. The band gap and the Rashba splitting near 
to the band edges of the pseudocubic CsPbG was found to be 
identical to those in (CH 3 NH 3 )Pbl 3 . Furthermore, the split¬ 
ting is fully suppressed in the centrosymmetric cubic CsPbG. 
This result suggests that precursors for the Rashba splitting 
are the strong spin-orbit interaction and distortions in the Pb- 
centred octahedron caused by the dipolar nature of CH 3 NH^ 
ions 10. A significant reduction of the Rashba splitting ob¬ 
served in (NH 4 )Pbl 3 perovskite structures El confirms va¬ 
lidity of this argument. 

The cubic structure is attractive for modelling due to its 
simplicity. However, we should keep in mind that the pseu¬ 
docubic structure is an idealization; the room temperature 
properties are governed by the tetragonal phase shown in 
Fig. [2ta). Distinctive features of the tetragonal structure are 
the presence of an octahedral tilting and CH 3 NH^ cation ro¬ 
tation. This leads to the question as to which extend the dis¬ 
cussion of the electronic structure of the pseudocubic phase is 
applicable to the tetragonal structure? In order to enable direct 
comparison between two electronic structures, the band struc¬ 
ture of low-symmetry tetragonal phase need to be mapped 
onto the pseudocubic Brillouin zone. This can be achieved by 
“unfolding” the supercell band structure. Since the tetragonal 
structure comprises of 2^ pseudocubic cells, each eigenstate 
in the supercell gives rise to eight k-points in the pseu¬ 
docubic Brillouin zone. For instance, r(0, 0, 0) unfolds into a 
set of eight k-points, three of which are obtained by permuta¬ 
tions of (1/2, 0, 0), another three additional k-points are per¬ 
mutations of (1/2,1/2,0), one k-point is at (1/2,1/2,1/2) 
and the T-point self. 

It is evident that T-point in tetragonal structure is cou¬ 
pled to R-point in the pseudocubic structure. However, it 
is not possible to say a priori what will be the strength of 
R-representation. Each unfolded k-point is assigned its own 
Bloch spectral weight w\^. The whole set satisfies normaliza¬ 
tion = 1- The spectral weight us used as an aug¬ 

mentation for the band structure plot F^(k, w\^) and allows to 
filter out spurious bands. The spectral weight is computed by 
projecting the supercell wavefunction onto a set of plane 
waves that are compatible with the wave vector k and recipro¬ 
cal lattice vectors of the pseudocubic Brillouin zone. Details 
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FIG. 2. Structure of /3-(CH3NH3)Pbl3 (a), /3-(CH3NH3)Pbl3 with a disorder in orientation of methylammonium molecules (b) and /3- 
(CH 3 NH 3 )Pb(Cla;Ii-cc )3 with X ^ 0.014 (c). The conventional tetragonal unit cell is outlined by a dashed contour on the panel (a). 


of this procedure are widely discussed in the literature ifT^ 

[BEMl. 

The effective band structure for tetragonal phase is shown 
in Fig.Oa). The band gap increases by approximately 0.1 eV 
when compared to the pseudocubic structure. The conduction 
and valence band edges of tetragonal phase retain 90 and 84% 
of the pseudocubic R-character, respectively. The F-character 
of the corresponding states is only 10% at most. The tetrago¬ 
nal phase exhibits a similar magnitude of the Rashba splitting 
as in the pseudocubic structure (Fig.fHb)). 

Now as the basis for comparison of the electronic struc¬ 
ture for high- and low-symmetry phases is established, we can 
proceed with disordered structures. A low energy barrier for 
rotation of the CHsNH^ ion and occurrence of multiple lo¬ 
cal minima |[^ leads to discretized reorientation dynamics 
of the C-N axis in (CH 3 NH 3 )Pbl 3 at the room temperature 
|[^ . The model for disordered orientation of methylammo¬ 
nium molecules was generated from /3-(CH3NH3)Pbl3 struc¬ 
ture (Fig. [2a)) by flipping the molecules 180° randomly with 
subsequent relaxation of atomic coordinates. The obtained 
structure is shown in Fig.|2b). It should be noted that the total 
energy of the disordered structure was even lower in compar¬ 
ison to the unperturbed tetragonal phase. 

The band structure of a disordered tetragonal phase is pre¬ 
sented in Fig. Ob). The disorder does not lead to signifi¬ 
cant changes in the band gap. The band edges retain their 
R-character at the level of 83 — 88%, which is almost identi¬ 
cal to that in the unperturbed tetragonal phase. Thus we can 
conclude that the orientational disorder of organic groups does 
not introduce localization effects at the band edges and does 
not hamper transport of charge carriers. This interpretation 
is based on the relation between the ambiguity of the Bloch 
character in reciprocal space and localization in the real space 
through the fundamental uncertainty principle. The only no¬ 
ticeable effect of the orientational thermal disorder is the re¬ 
duced Rashba splitting, which can make the spin selection 
rule, which promotes dissociation of optical excitations, less 
efficient. 

Finally, we turn the discussion to a compositional disorder 
in mixed halide perovskite structures. Even though the Cl in¬ 


corporation in an iodide-based structure is relatively low (be¬ 
low 3-4%) lH^, the recombination time and diffusion length 
of optical excitations are prolonged by 1-2 orders of mag¬ 
nitude The origin for a superior performance of 

(CH 3 NH 3 )Pb(Cla;Ii_a ^)3 perovskite thin film solar cells re¬ 
mains an open question. The plausible scenarios include 
Cl“ ions involved in the following functions: facilitation the 
release of excess CH 3 NH^ during annealing (H], passiva¬ 
tion of grain boundaries and the interface with the substrate 
in, 113], creation of an in-homog^eous built-in electric field 
profile similar to p-i-n junction and reduction of the 

carrier effective mass in structures with high Cl content lld^ . 

In order to model a dilute substitutional 
(CH 3 NH 3 )Pb(Cla;Ii_a ^)3 alloy, one of 24 I-atoms in the 
supercell was replaced by Cl that resulted in the effective 
composition of x ~ 0.014. The structural relaxation in the 
vicinity of chlorine reveals asymmetric Pb-Cl bonds di < d 2 
(see Fig. I2c)). The Pb-Cl bond lengths of 2.75 and 3.48 A 
deviate significantly from the Pb-I equilibrium bond length of 
3.25 A. The shorter Pb-Cl bond length agrees with the sum of 
Pb^+ and Cl“ ionic radii, 0.94 and 1.81 A respectively lld^] . 
The asymmetric configuration could be an efficient way to 
mediate local strain due to the large size mismatch between 
cations and can be a driving force for segregation of Cl at 
interfaces. 

The effective band structure of 

(CH3NH3)Pb(Clo.oi4lo.986)3 IS shown in Fig. Oc). In- 
corporation of chlorine in dilute quantities does not change 
the band gap, which is in accord with observations made 
by Colella et al. HH]. The Rashba splitting and the Bloch 
character of the valence and conduction band edges remain 
intact with the undopped (CH 3 NH 3 )Pbl 3 . Perturbations due 
to Cl are observes for states within the conduction and valence 
bands as pointed out by arrows in Fig.Oc). Since the access 
energy required to reach perturbed states is greater than the 
thermal energy, those perturbations in the band structure will 
unlikely have an adverse effect on transport coefficients of 
the mixed halide perovskites. 
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FIG. 3. Effective relativistic band structure of /3-(CH3NH3)Pbl3 (a), /3-(CH3NH3)Pbl3 with a disorder in orientation of methylammonium 
molecules (b) and /3-(CH3NH3)Pbl3(i_aj)Cl3cc with x ^ 0.014 (c). The size and colour of the markers correspond to the spectral weight of 
supercell eigenstates projected to the Bloch states in the pseudocubic Brillouin zone (see inset in Fig.lHb)). 


IV. CONCLUSIONS 


The relativistic band structure of low-symmetry p- 
(CHsNHsjPbls tetragonal phase is mapped onto the pseu¬ 
docubic Brillouin zone using an unfolding technique. In spite 
of significant differences in the structure (the presence of an 
octahedral tilting and CHsNH^ cation rotation), the band 
structures of both phases share many similarities. A strong 
coupling is observed between F-point in tetragonal structure 
and R-point in the pseudocubic structure. The conduction and 
valence band edges of tetragonal phase retain 80-90% of the 
pseudocubic R-character. The V -character of the correspond¬ 
ing states is very weak in contrast to the common expectation. 
Both pseudocubic and tetragonal structures exhibits a simi¬ 
lar magnitude of the Rashba splitting, which plays an impor¬ 
tant role in the long carrier lifetime due to constrains imposed 
by the spin section rules. At room temperature, the Rashba 
splitting can be reduced due to an orientational disorder of 
CHsNH^ cations. At the same time, the disorder neither af¬ 
fects the band gap nor reduces the R-character of the band 


edges and thus does not hamper transport of charge carriers. 

Mixed halide (CH 3 NH 3 )Pb(Cla^Ii_a ^)3 structures with di¬ 
lute concentration of chlorine reveal existence of asymmetric 
Pb-Cl bonds as a way to mediate local strain due to the large 
size mismatch between cations that can also be a reason for 
segregation of Cl at interfaces. The band gap, Rashba splitting 
and the Bloch character of the valence and conduction band 
edges in perovskite structures with dilute substitutional chlo¬ 
rine remain intact with the undopped (CH 3 NH 3 )Pbl 3 . Pertur¬ 
bations due to the incorporated chlorine are observes for states 
located energetically far from the band edges. This rules out 
the possibility of adverse effect of chlorine substitutional in¬ 
corporation on transport coefficients of the mixed halide per- 
ovskites. 
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